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ABSTRACT 


A  new  type  of  moments  has  been  achieved  by  substituting  in 
the  central  moments  the  smallest  value  of  the  sample  for  its  mean. 

The  new  moments  have  the  same  advantage  as  the  central  moments  of 
being  independent  of  the  location  parameter  but  for  certain  values 
of  the  shape  parameter  they  have  less  variance  and  thus  are  pre¬ 
ferable  for  estimating  purposes.  The  asymptotic  properties  of  four 
estimators,  three  of  them  composed  of  the  new  moments  and  one  of 
them  of  central  moments  have  been  examined.  It  could  be  concluded 
that  for  the  shape  parameter  cLi  0. 5  the  estimator,  which  was  composed 
of  the  first  and  second  order  moments  of  the  new  type,  was  by  far  the 
most  efficient  one.  Small-sample  properties  of  the  new-moments 
estimators  have  been  appraised  by  use  of  extensive  Monte-Carlo 
studies  and  it  could  be  stated  that  the  same  conclusion  applies  also 
to  small  and  moderate  sanple  sizes. 


Distribution  of  this  abstract  is  unlimited.  It  may  be  released 
to  the  Clearinghouse,  Department  of  Commerce,  for  sale  to  the  general 
public. 
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Moments  about  Smallest  Sample  Value 


1.  Introduction 


The  central  moments,  defined  by 


my-E  (xi  -  x)v/« 


where  T  denotes  the  mean  of  the  sample  values,  has  the 
useful  property  of  being  independent  of  the  location  para¬ 
meter  as  easily  seen  by  introducing  the  standardized  variate 

z  -  (x-p)/p 


x  ■  P  z  +  n 


pV£(  t±  -  S)7n-  pV  .  fv(«)/n 


is  uniquely  determined  by  a  and  p. 

Considering  that  the  variances  of  the  mean  x  and  of 
the  smallest  sample  value  x^  from  a  Weibull  population  are 

2a 

Var (f)«ji2/n  and  Varfx^)  - fig/n 

it  can  be  oonoluded  that_for  a  >0.5  the  statistio  x^  is 
much  more  precise  than  x  for  large  sample  sizes  n. 

Thus  it  seems  plausible  that  the  efficiency  of  the  moment 
estimators  oan  be  increased  by  substituting  x^  for  If  in 
equ.(l). 

In  this  way  a  new  type  of  moments  is  obtained,  which  will 
be  called  "moments  about  smallest  sample  value",  be  denoted  by 
r  and  defined  by 


rv  *  £(xi  ~  xi)Vn 

The  properties  of  these  moments  and  their  usofulneas  for 
estimating  purposes  will  now  be  examined. 

2.  Properties  of  the  new  moments 

Introducing  (2)  into  (5)  we  have 

rr  -  pTE(si  -  *1)T/«  -  PT  f^  («)/n 
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from  which  it  can  be  concluded,  since  z  is  uniquely  deter¬ 
mined  by  a,  that  the  moments  r  ,  just  like  the  moments  m  , 
are  independent  of  the  looation  parameter  p.  r 

It  should  be  noted  that,  by  definition,  the  first  oentral 
moment 

m1  «  l(xi  -S)/n  «  0 

whereas  the  moment 

r!  -  Zix^xj/n  -  x  -  ^  >0 
and  thus  offers  a  new  tool  for  estimating  purposes. 

Prom  (8)  it  follows  that  the  expected  value  of  r^ 

Er^-Ex-  E 

and  considering  that 

B  x  ■  +  p  g^  and  E  x^  ■  p  +  p  gj/n“ 

where,  using  a  notation  proposed  by  Dubey, 

gj^  -  «  1  and  g^,-  (va)! 
the  expected  value  of  r^  becomes 

Erx  -  pgx(l  -  n~“) 

Poxther,  from  equ.{8) 

Varfr^)  •  Var(x)  ♦  Var(x^)  +  2 Cov  (J,!^ 

The  first  two  terms  are  given  by  equ.(4)t  while  the 
third  term  inoludes  terms  of  the  type  E(x^x  ),  whioh  can 
be  computed  by  means  of  formulae  derived  by  ^  Lieblein  [l]. 

In  the  same  way,  the  expeoted  value  of  Er^  of  the 
eeoond  moment 

2  2 
r2"r^i  ”  Xl^  &2~  2x\'*l  *  *1 

can  be  evaluated. 

The  expeoted  values  and  the  variances  of  the  moments  r^ 


take  Tory  simple  forms  for  large  sample  sizes,  since 

ry  ■  l(xi  -  rL)Vn-  pVl(zi  -  - #■  pV£z’'/n  (15) 

as  n  —  >  ®  . 

From  formulae  given  in  Scientific  Report  No. 6  (SR  6) 
it  then  follows  that 

Erv-pV.gv  (16) 

and 

n.Var(ry)  -  p2v(g2y- g^2)  and  n.Cov(rv,rw)  -  pV+W(gv+w~  (*7) 

3.  Use  of  the  moments  for  estimating  shape  parameters 

From  equ.(6)  it  is  easy  to  see  that  by  use  of  two  moments 
r  of  different  orders,  an  expression  can  be  constructed  that 
i?  independent  of  p  and  thus  a  function  of  a  only,  for  ex¬ 
ample,  the  expressions 

fi12"r2/rl  ®13 "  r3^rl  H23"r32^r23  (l8) 

Each  such  expression  oan  be  used  as  an  estimator  of  a. 

Since  a.  -0,  the  only  corresponding  estimator,  oooposed 
of  oentral  moments,  is 

M23  "  {19) 

The  efficiency  of  these  estimators  will  now  be  examined. 

Due  to  the  fact  that  the  moments  r  for  large  sample  sizes 
tend  to  the  moments  about  origin  oithe  standardized  variate, 
closed  expressions  oan  in  this  case  be  deduced.  For  small  and 
moderate  sample  sizes  the  Monte-Carlo  Method  has  to  be  used. 

4.  Asymptotio  properties  of  various  estimators 

If  n  — — w  Co,  the  estimators  R  tend  to 

*12-*A2  B23-*32/*23  <*-0)  (20) 

where  a^  are  the  moments  about  origin  of  the  standardised 
variate. 
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According  to  a  theorem  stated,  and  proved  by  CraaAr  [2,p.352] 
for  the  case  of  a  function  H(m  ,m  )  of  two  oentral  moments, 
we  have,  denoting  by  H  ,  H-  aXd  w  the  values  assumed  by 
the  function  B  and  it§  first  order  partial  derivatives  in  the 
point  m  *fi  ,  a  ■  fi  ,  the  mean  and  the  variance  of  E  are 
given  byv  v  w 

E(H) -Ho  +  0(l/n)  (21) 

Var(H)  -  H*.  Var(r?v)  +  2H1H2  .  Cov(«V“w)  +  H2  *  *“(%)  +  0(l/n3/2) 

Since  these  formulae  are  valid  for  any  mean  a  of  the  po¬ 
pulation,  including  m-0,  they  are  valid  also  for  the  moments 


4*1  Mean  and  variance  of  the  estimator  _ 

Prom  equ.(20)  we  thus  have 

H^a^-a^  (|x -0) 

and 

'  Hl“'2«2/Pgl3  '  V1/**2*! 

and  by  equ.(l?) 

n.Var(r^)  -  p2(«2  -  g^)  |  n.  Cov^,^)  -  p3(g3  -  g1g2)  , 

n  .  Var(r2)  -  P4(g4-g2  ) 
which  after  some  calculations  results  in 

B(R12)  -g^/gj2 

aiid 

n.VartR^)  -  (4gg  +  -  K 

If  now  a  value  R-2  has  been  oomputed  from  a  random 
sample,  the  estimated  value  a*  of  the  parameter  a  is 
obtained  by  equating  the  actual*  value  R.-  to  the  expec¬ 
ted  value  ECR^). 

Hence 

H12  -  g^gj2  •  f(a*2) 


(22) 

(23) 

(24) 


(25) 

(26) 


(27) 
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Some  values  of  the  function 


f(a)  are  presented  in  Table  1. 


Since  VarfR^g)  — 0  as  n  - *>  ®  and  R^2  - -  f(«o) 

where  a  is  the  true  value  of  the  parameter,  it  follows  that 
the  valuS  a.  2  computed  from  equ«(27)  is  a  consistent  estimate 
of  a.  PUrtner,  that  the  tangent  curve  at  a  can  be  substi¬ 
tuted  for  the  function  f(a^2).  0 


Bence 


f(«0)  ♦  («*2  -  «0) .  f><«0)  .  S12 

(28) 

and 

Var(o£2)  -  Var(H12)/[f»(ao)]2 

(29) 

From 

equ. (14)  of  SR  7 

dg^/da-  v  .  g^  X"(va) 

(30) 

where 

X"  (a)  -  the  digamma  function,  which  is  tabulated. 

and 

^(a)  •d(g2/g2)/da  -  2gg[  (2«e)  -  X'(a)]/g12 

(31) 

The  asymptotic  efficiency  thus  becomes 
n.AsyEff (a*2)  -  [gg/g2  +  g^Ag*  -  "  1/^/C  f  (2o)  -/'(a)]2  (32) 

for  a  >  0.5. 

4»2  Means  and  variances  of  the  estimators  and  Rg^ 

Proceeding  in  the  same  way  as  above  the  asymptotic  effici¬ 
encies  are  found  to  be 

n.AayEfffrjjHgg/g^  +  g6/9g2  -  4 /9]/[  f  ($«)-/' (  a)]2  (33) 

and 

n.AsyEff(«g3)-[g494«2+  g^g2  -  g^/ig^  -  l/36]/[  r  (3o)-r(2a)]2  (34) 
Values  of  these  variances  are  presented  in  Table  2. 
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which  values  are  presented  in  Table  3  of  SH7  for 
v  ■  2*6,  we  have 

E(M23)  -  k3/k23/2 

and 

n.Var(M23)-k6/k^+  2.25k2k4A25  +  8.75^/^+  9.0  + 

-  3k3k5/k24  -  6k4/k22 

Some  values  of  E(M_,)  are  presented  in  Table  1,  and 
of  Var(ll03)  in  Table  2. 


4.4  Efficiencies  of  the  estimators 

The  lower  limit  attainable  of  the  variances  of  the  esti¬ 
mate  of  a,  denoted  by  Vara,  is  given  in  SR  4  for  various 
alternatives.  Some  values  for  the  case  that  all  three  para¬ 
meters  are  unknown  are  presented  in  Table  2. 


(35) 

(36) 

(3?) 


The  efficiency,  as  defined  by  R.A.  Fisher,  of  any  estimate 
is  obtained  by  dividing  Var  o  by  the  variance  of  the  estimate 
in  question.  Some  values  resulting  from  this  procedure  are 
presented  in  Table  3,  from  which  is  seen  that  all  the  estima¬ 
tors  R  are  more  efficient  than  M?,  and  that  R--  is  the 
most  efficient  of  them  all.  It  should  be  noted  that  this 
conclusion  is  valid  for  a  >0.5  only.  For  a  <0.5  the 
estimator  M„3  is  more  efficient  then  the  estimator  R. 
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5.  Small-sample  properties  of  the  estimators  and 

5.1  Monte-Carlo  study  of  the  distributions  of  the  estimators 

An  extensive  Monte-Carlo  study  has  been  performed  in  order  to 
determine  the  sampling  distributions  of  the  two  estimators.  By  use 
of  an  IBM  7090  computer  10,000  random  samples  for  each  combination 
of  o  -0.001,  0.01,  0.1,  0.3,  0.5,  0.7,  0.9,  1.0,  1.5,  2,0,  5.0 
and  10.0  and  n-3,  4,  5,  10  and  20  have  been  produced,  f^om  which 
the  values  of  B^p  and  Mp,  have  been  computed  and  classified  in 
40  classes.  The1  means  and  the  standard  deviations  of  the  ten 
thousand  values  of  the  sixty  combinations  have  at  the  same  time  been 
determined.  The  values  of  are  presented  in  Tables  4,  5  and  6 

for  n-5,  10  and  20,  and  those  of  M_,  in  Tables  7,  8  and  ')  for 
n  -  5,  10  and  20.  J 

Since  the  moments  m  and  r  tend  to  zero  as  a  tends  to 

zero,  the  estimators  tend"to  the  ratio  0/0.  For  th?s  reason, 

the  data  corresponding  to  a  -0.001  have  been  omitted  as  being 
somewhat  unreliable,  and  so  Slso  the  data  for  a  -5-0  and  10.0 
as  being  of  less  practical  interest.  0 

The  cumulative  distributions  of  and  have  been 

plotted  as  demonstrated  in  Figs.l  and  2  , ^respectively. 


5.2  Means,  medians  and  modes  of  the  estimators 

The  Monte-Carlo  determined  means  B1?  and  M?,  ,  computed  from 
the  10,000  values,  are  presented  in  Tables  10  and  11  ,  including 
also  the  values  for  oo»5*0  and  10.0. 

It  was  found  that  the  statistics  R.p  ,  R. ^  and  R^p  and  corre¬ 
sponding  statistics  with  regard  to  Mp,  could  with  excellent  appro¬ 
ximation  be  expressed  as  polynomials  ^of  the  second  degree  in  a  . 
Fitted  values  and  corresponding  equations  are  presented  in  Tabled 
12  and  13. 

The  medians  R^p  and  the  modes  R^p  have  been  graphically 
evaluated  from  the  distribution  curve!  in  Fig.l  and  the  corre¬ 
sponding  curves  for  n-5  and  n-10.  These  values  are  also  in¬ 
cluded  in  Table  12. 
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The  Monte-Carlo  deterrined  variances  n.Var(R,2)  and  n.Var^ 
obtained  by  squaring  the  standard  deviations  and  multiplying  the 
squares  by  n,  are  presented  in  Tables  14  and  15* 


Variances  of  the  estimators 


Varianoes  of  the  estimates 


On  the  condition  that  the  estimate  at ,,  is  obtained  by  equa¬ 
ting  observed  values  of  R_2  to  the  §xp|cted  value  of  R12,  as 
indicated  by  equ.(27),  the  variance  of  a.pis  computed  from  the 
variance  of  R12  by  use  of  equ.(29).  This  procedure  is  strictly 
oorrect  for  n-oo  ,  but  an. error  will  be  introduced  for  finite 
sample  sizes*  Nevertheless,  this  method  will  be  applied,  and  the 
result  will  be  used  as  a  measure  of  the  variance. 


From  the  equations  in  Tables  12  and  13  it  follows  that  for 


n  «  5  dl12 /da -0.2835  +  0*4280  a 
n  -  10  dl12 /da- 0.2321  +  0.9324a 
n-20  d?12 /da -0.1711  +  1.2724a 


n  -  5  dM2;} /da -1.050  -  0.480  a 

n  -  10  dM23  /da  -  2.094  -  1.152  a 
n-20  dfi23 /da- 2.830  -  1*552  a 


+  By  us^of  these  formulae  and  equ.(29)  the  variances  of 
a.p  and  a2  have  been  computed.  The  results  are  presented 
inTable  lb.  ^t  is  of  interest  to  note  that  for+a>0.5  the 
variance  of  a23  is  always  larger  than  that  of  a^2  . 


It  is  obvious  that  this  method  of  estimation  will  prov^Le 
estimates  with  negative  bias  of  a^2  and  positive  bias  of  a* 
but  these  biases  will  tend  to  zero  with  n  .  However,  if  a 
large  set  of  samples  are  used  for  computing  and  the  mean 

of  the  values  R^2  is  introduced  into  eq  .(27),  an  unbiased 
estimate  will  result. 


The  above-mentioned  method  of  estimation  is  not  the  only 
one  possible.  It  would  be  quite  reasonable  to  equate  the  actual 
value  of  fi12  to  the  median  S^2  .  The  advantage  of  doing  so  would 


be  that  a  correct  median  estimate  will  be  obtained,  that  is, there 
would  be  50  $  probability  of  having  a  value  larger  than  and  50$ 
probability  of  having  a  value  smaller  than  the  true  value  of  a. 
Most  likely,  the  variance  of  this  estimate  will  be  a  little  lar¬ 
ger,  but  the  bias  a  little  smaller  than  those  of  the  preceding 
method. 


A  third  alternative  of  estimation  would  be  to  equate  the 
actual  value  of  IL  to  the  mode  R.  p  •  Such  an  estimate  could 
be  called  the  most  probable  estimate!  To  illustrate  this  state¬ 
ment  by  an  example,  suppose  that  from  a  sample  of  size  n *  20  a 
value  R12«1.45  has  been  computed.  Knowing  nothing  about  the 
true  value  of  a,  it  may  be  assumed  that  there  is  the  same  pro¬ 
bability  for  all  values  of  a  between  0  and  1  .  (In  many  cases 
these  will,  for  logical  reasons,  be  the  limits).  Then  it  will 
be  found  from  Table  6  that,  the  total  number  of  R.?  being  7422, 
the  probabilities  of  a  value  I.45  coming  from  a  population 
with  the  shape  parameter  a  will  be  for 


a  -  0.1,  0.3,  0.5,  0.7,  0.9 
P  -  5.7  16.8  37.1  25.6  7.9  $ 


The  confidence  limits  of  this  estimate  are  easily  set. 

Very  likely,  the  variance  of  this  estimate  will  be  somewhat 
larger  than  those  of  the  two  preceding  methods,  but  possibly 
the  bias  will  be  the  least  one. 

It  should  be  noted  that  the  differences  between  the  three 
methods  will  disappear  with  increasing  sample  size,  since  all 
the  distributions  are  tending  to  the  same  normal  distribution. 

5.5  Efficiency  of  the  estimates 

In  order  to  obtain  the  efficiency  of  the  various  estimates, 
the  varianoe  Var 0  has  to  be  divided  by  the  respective  variances. 
The  result  thus  reached  is  presented  in  Table  17.  It  should  be 
observed  that  the  indicated  values  of  the  lower  limit  n.Varfi 
is  strictly  valid  for  large  sample  sizes  only,  so  another  appro¬ 
ximation  has  thereby  been  introduced.  However,  it  is  believed 
that  the  values  will  be  sufficiently  accurate  for  a  comparison 
between  the  estimates,  even  if  they  do  not  give  the  real  effi¬ 
ciencies,  whioh  may  be  better  than  indicated  in  the  table. 
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Conclusions 


An  examination  of  the  asymptotic  properties  of  four  esti¬ 
mators,  three  of  them  based  on  moments  about  the  smallest  sample 
value  and  one  on  central  moments,  has  proved  that  for  a  >0.5 
the  first  type  of  estimators  is  much  more  efficient  than  the 
second  type,  and  that  the  estimator  composed  of  the  first  and 
the  second  moment  about  smallest  sample  value  is  by  far  the  most 
efficient  one. 

Encouraged  by  this  result,  an  extensive  Monte-Carlo  study 
has  been  performed  with  the  result  that  the  same  conclusion  ap¬ 
plies  also  to  moderate  and  small  sample  sizes. 

The  computational  work,  when  using  any  one  of  the  moment 
estimators  is  considerably  less  than  that  required  by  the  Maxi¬ 
mum  Likelihood  Method,  which,  however,  has  the  advantage  of  being 
more  efficient. 
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Table  1.  Expected  values  of 


and  M„,  for  n  ■  oo 


ft 

e(m23) 

a 

E<V 

E<V 

0.00 

1.0000 

-1.1396 

0.60 

1.3801 

0.8960 

0.01 

1.0002 

-1.0809 

0.70 

1.5045 

1.1604 

0.10 

1.0145 

-0.6376 

0.80 

1.6480 

1.4295 

0.20 

1.0524 

-0.2541 

0.90 

1.8124 

1.7080 

0.30 

1.1093 

O.O687 

2.0000 

2.0000 

0.40 

1.1831 

0.3586 

3.3953 

3.8023 

0.50 

1.2732 

0.6311 

6.0000 

6.6188 

Table  2.  Asymptotic  variances  of  various  estimates 


n 

n.Var 

n.Var  ^ 

n.Var 

n.Var  «*3 

n.Var 0 

0.5 

0.15574 

0.16715 

0.20545 

0.88939 

0.15198 

0.6 

0.23484 

0.27000 

0.36939 

1.37268 

0.21865 

0.7 

0.34454 

0.42751 

0.65061 

2.04212 

0.29789 

0.8 

0.49676 

0.66876 

1.12759 

3.33016 

0.38908 

mm 

0.70731 

1.03759 

1.92656 

5.17201 

0.49242 

Kb 

1.00000 

1.60000 

3.25000 

8.00000 

0.60793 

Table  3.  Asymptotic  efficiency  of  various  estimates 


a 

f  °1 2 

Eff 

Eff  o^3 

Eff«»3 

0.5 

97-59 

90.92 

73-97 

17.09 

0.6 

93.19 

81.06 

59.25 

15.94 

0.7 

86.46 

69.68 

45-78 

14-59 

0.8 

78.32 

58.18 

34.50 

11.68 

0.9 

69.62 

47.46 

25.56 

9.52 

1.0 

60.79 

38.00 

18.71 

7.60 
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Table 


4.  Claes  frequencies,  means  and  star  dard  deviations  of  ftj... 
for  aample  aise  a  ■  5  _ 


Tipper 

olass 

limit 

The  true  shape  parameter  a 

0 

equal  to 

.01 

.1 

.3 

•5 

.7 

•9 

1.0 

1.5 

2.0 

in 

0 

0 

0 

0 

0 

0 

0 

0 

0 

mwvm 

1616 

1263 

728 

386 

228 

157 

131 

55 

28 

mm 

2641 

2493 

1891 

1327 

907 

650 

529 

241 

138 

1.5 

1735 

1833 

1780 

1.554 

1134_ 

815 

729 

382 

202 

1.6 

1226 

1304 

1535 

1418 

1319 

1028 

910 

467 

277 

1.7 

980 

1028 

1185 

1329 

1274 

1148 

1051 

697 

442 

1.8 

609 

669 

842 

1041 

1063 

1050 

1018 

675 

498 

1.9 

369 

424 

566 

740 

882 

876 

824 

695 

468 

2.0 

252 

303 

417 

550 

685 

784 

811 

687 

521 

2.1 

1§4 

209 

279 

409 

547 

652 

<j&T 

“644 

500 

2.2 

136 

145 

239 

295 

428 

518 

582 

620 

533 

2.3 

74 

110 

146 

269 

348 

429 

476 

561 

523 

2.4 

55 

61 

114 

171 

258 

352 

374 

558 

543 

2.5 

33 

45 

75 

127 

203 

308 

292 

472 

532 

2.6 

27 

28 

54 

78 

153 

210 

297 

459 

2.7 

14 

22 

36 

67 

120 

185 

204 

349 

362 

2.8 

11 

16 

25 

73 

95 

136 

180 

285 

326 

2.9 

7 

8 

23 

34 

63 

142 

144 

225 

290 

3.0 

5 

9 

12 

34 

62 

99 

129 

210 

289 

3.1 

7 

6 

9 

17 

60 

93 

121 

202 

231 

3.2 

4 

6 

17 

22 

42 

62 

95 

178 

257 

3.3 

3 

4 

1 

12 

28 

62 

69 

168 

229 

3.4 

5 

4 

8 

13 

23 

55 

70 

155 

207 

3.5 

2 

4 

4 

9 

21 

57 

156 

188 

3.6 

2 

2 

7 

5 

12 

37 

63 

124 

173 

3.7 

1 

1 

2 

6 

17 

24 

38 

123 

195 

3.8 

0 

1 

2 

5 

7 

17 

31 

114 

167 

3.9 

0 

0 

0 

4 

7 

21 

28 

97 

156 

4.0 

0 

0 

1 

2 

1 

15 

19 

108 

145 

® 

2 

2 

2 

13 

2b 

48 

384 

~TT?T 

H12 

1-5174  1 

.5443  1 

.6166 

1.7075  1.8163  1-9395 

2.0051 

2.3469  2.6718 

S.D. 

.26665  • 

28020  . 

31659 

.36469 

.42678 

.50075 

.54039 

.73342 

.87908 

12 


The  true  ahape  parameter  a  equal  to 


.3  .5 


H1?  1.3108  1.3399  1.4245  1.5426  1.6993  1.8959  2.0083  2.6752  3-4156 
3.D.  .14558  .15588  .18590  .23389  .31210  .42749  49847  .93100  1. 360 31 
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Table  6.  Class  frequencies,  means  and  standard  deviations  of  -R12 
for  sample  si  as  n  -  20 


Upper 

The  true 

shape  parameter  «o 

equal  to 

olass 

limit 

.01 

.1 

.3 

.5 

.7 

.9 

1.0 

1.5 

m 

1.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1.1 

799 

5 

0 

0 

0 

0 

0 

0 

1.2 

4947 

4191 

1297 

109 

3 

0 

0 

0 

0 

1.3 

3009 

3624 

3857 

1302 

I65 

15 

6 

0 

0 

1.4 

947 

1356 

2938 

2993 

921 

166 

69 

1 

0 

1.5 

240 

420 

1245 

2753 

1900 

588 

276 

_ 7 _ 

0 

1.6 

45 

93 

454 

T5&f 

2302 

1137 

S84 

25 

3 

1.7 

11 

24 

147 

733 

1806 

1534 

1053 

90 

4 

1.8 

1 

4 

38 

329 

1244 

1535 

1274 

156 

12 

1.9 

0 

1 

14 

138 

701 

1370 

1320 

259 

43 

2.0 

1 

1 

2 

42 

398 

1058 

1191 

421 

48 

2.1 

0 

0 

1 

24 

283 

793 

955 

479 

80 

2.2 

0 

0 

1 

7 

124 

534 

799 

528 

142 

2.3 

0 

0 

1 

4 

77 

376 

582 

622 

172 

2.4 

0 

0 

0 

0 

22 

268 

453 

633 

233 

2*5 

0 

0 

0 

1 

17 

194 

324 

659 

265 

2.6 

0 

0 

0 

2 

7 

135 

~23B~ 

561 

294 

2.7 

0 

0 

0 

1 

12 

96 

205 

581 

290 

2.8 

0 

0 

0 

0 

5 

59 

136 

547 

354 

2.9 

0 

0 

c 

0 

2 

49 

100 

460 

365 

3.0 

0 

0 

0 

0 

3 

27 

90 

402 

360 

3.1 

0 

0 

0 

0 

1 

17 

a 

399 

375 

3.2 

0 

0 

0 

0 

0 

8 

45 

332 

373 

3.3 

0 

0 

0 

0 

1 

9 

26 

375 

3.4 

|  0 

0 

0 

0 

0 

8 

25 

362 

3.5  0 

0 

0 

0 

0 

2 

.  19 

238 

332  . 

3.? 

i  0 

0 

0 

0 

0 

4 

16 

210 

289 

3.7 

!  o 

0 

0 

0 

0 

5 

6 

193 

315 

3.8 

!  0 

0 

0 

0 

0 

3 

4 

178 

294 

3.9 

i  0 

0 

0 

0 

1 

0 

7 

121 

268 

.  4.0 

0 

0 

0 

0 

0 

4 

6 

142 

241 

ID 

!  0 

j 

0 

0 

0 

0 

T“ 

25 

liW 

41  oF 

j 

\2  |1.1996 

1.2274 

1.3133 

1.4411 

1.6211 

1.8605 

2.0038 

2.9448 

4.1628 

!  S.D. 
1 _ 

,.08691 

1 - 

•  09347 

.11207 

.14486 

•  209C  8 

.32129 

.39863 

.99530 

1.80312! 
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Table  7.  Class  frequencies,  means  and  standard  devlaticne  of 
for  sample  size  n  «  5 


Upper 

The  true  shape  parameter  «o 

equal  to 

0l&8B 

limit 

.01 

.1 

.3 

•  5 

.7 

.9 

1.0 

1.5 

2.0 

am 

16 

0 

0 

0 

0 

0 

0 

0 

0 

-  1.8 

1 

0 

0 

0 

0 

0 

0 

0 

0 

-  1.6 

5 

0 

0 

0 

0 

0 

0 

0 

0 

-  1.4 

124 

85 

44 

20 

15 

11 

10 

4 

2 

-  1.2 

525 

411 

179 

108 

61 

34 

26 

14 

10 

-  1.0 

615 

550 

348 

!?6 

100 

7  2 

58 

25 

11 

ETK-1 

740 

639 

392 

286 

176 

114 

96 

33 

21 

837 

830 

595 

366 

269 

182 

150 

78 

31 

-  0.4 

998 

1020 

875 

637 

453 

320 

264 

148 

82 

-  0.2 

1325 

1316 

1199 

986 

780 

638 

595 

326 

234 

0.0 

1303 

1306 

1265 

1104 

8G9 

723 

639 

412 

263 

MEM 

911 

1105 

1273 

1292 

1141 

944 

867 

570 

384 

BO 

891 

1015 

1226 

1444 

1479 

1397 

1349 

1034 

739 

WM 

557 

642 

945 

1066 

1219 

1309 

1277 

1260 

1260 

0.8 

430 

422 

596 

812 

973 

1044 

1079 

923 

852 

1.0 

-263__ 

306 

445 

677 

789 

?1} 

978 

1031 

842 

1.2 

197 

206 

346 

516 

763 

694 

929 

1200 

1105 

1.4 

127 

121 

213 

404 

683 

985 

1110 

1447 

1649 

1.6 

56 

26 

59 

104 

210 

420 

573 

1495 

2515 

1.8 

22 

0 

0 

0 

0 

0 

0 

0 

0 

2.0 

10 

0 

0 

0 

0 

0 

0 

0 

0 

® 

47 

0 

0 

0 

0 

0 

0 

0 

0 

-.1948 

-.1732 

0148 

.1924 

.3486 

.4812 

•  5392 

.7630 

.9099 

S.D. 

1.07705  .61352 

.60428  .60034 

.59812 

.59401 

.59089 

.567  18 

.53784 
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Table  8.  Class  frequencies,  means  and  standard  deviations  of 


for  sample 

sise  n » 

JLO 

0— 

Upper 

class 

limit 

The  true  shape  parameter  ctQ 

equal  to 

.01 

.1 

.3 

.5 

.7 

.9 

1.0 

1.5 

2.0 

-  2.0 

123 

34 

3 

0 

0 

0 

0 

0 

0 

-  1.8 

131 

75 

4 

1 

0 

0 

0 

0 

0 

-  1.6 

235 

106 

9 

2 

1 

0 

0 

0 

0 

-  1.4 

334 

190 

37 

6 

2 

1 

0 

0 

0 

-  3.2 

514 

320 

71 

10 

2 

1 

1 

0 

0 

i-  1.0 

mkim 

561 

165 

32 

10 

5 

2 

0 

0 

947 

832 

303 

79 

20 

5 

S~ 

0 

0 

aiwi 

1090 

1055 

554 

201 

68 

28 

13 

5 

0 

—  0.4 

1278 

1249 

890 

386 

161 

62 

46 

9 

0 

-  0.2 

1288 

1378 

1251 

691 

321 

162 

113 

25 

13 

0.0 

1201 

1292 

1503 

1105 

617 

329 

246 

60 

19 

0.2 

8l9 

1111 

1466 

1399 

1002 

637 

500 

158" 

60 

0.4 

573 

791 

1295 

1501 

1259 

927 

796 

335 

139 

0.6 

342 

487 

997 

1354 

1434 

1219 

1053 

553 

314 

0,3 

197 

288 

676 

1176 

1348 

1342 

1269 

857 

498 

1.0 

81 

115 

386 

853 

1185 

1284 

1294 

1008 

774 

1.2 

35 

62 

185 

503 

912 

1151 

1221 

1689 

877 

1.4 

23 

33 

112 

316 

647 

916 

1011 

1226 

1057 

1.6 

10 

10 

47 

I84 

391 

665 

767 

1104 

1167 

1,8 

2 

8 

28 

117 

278 

461 

571 

789 

938 

—2,0. 

0 

2 

11 

46 

177 

341 

mtm 

786 

809 

2.2 

2 

0 

6 

28 

107 

255 

330 

696 

'  851 

2.4 

0 

1 

0 

8 

48 

148 

218 

625 

935 

2.6 

0 

0 

1 

2 

9 

58 

102 

56I 

1079 

2.8 

0 

0 

0 

0 

1 

3 

8 

114 

00 

0 

0 

0 

0 

0 

0 

0 

0 

0 

m 

-.4976 

-.3324 

.0363 

.3670 

,6478 

.8821 

.9840  1.3768 

1.6356 

1 

.63159 

.59231 

.54587 

.54962 

•57581 

.60266 

.61371 

.64178 

.63899 
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Table  9.  Claaa  frequencies,  means  and  standard  devi&tiona  of  ML 
for  sample  size  n  »  20 


The  true  shape  parameter  a  equal  to 


.3  .5 


5 

0 

34 

3 

113 

25 

254 

73 

-.7315  -.4630  .0471  .4830  .8569 

1.1807  1.3263 

1.9187 

2.3368 

.57480  .49993  .41990  .45010  .52099 

.59586  .63010 

.75271 

.80651 

n«  3 

n*4 

n  *  5 

n  *  10 

n-  20 

1*7756 

1.6165 

1.5174 

1.3108 

1.1996 

1.7919 

1.6402 

1.5443 

1.3399 

1.2274 

1.8320 

1.7022 

1.6166 

1.4245 

1.3133 

1.8773 

1.7766 

1.7075 

1.5426 

1.4411 

1.9261 

1.8617 

1.8163 

1.6993 

1.6211 

1.9767 

1.9545 

1.9395 

1.8959 

1.8605 

2.0021 

2.0026 

2.0051 

2.0083 

2.0038 

2.1238 

2.2441 

2.3469 

2.67  52 

2.9448 

2.2299 

2.4643 

2.6718 

3.4156 

4.1628 

2.5807 

3.2065 

3.8012 

6.4644 

10.6992 

2.7700 

3. 5861 

4.3714 

8.1353 

14.9989 

Table  11.  Monte-Carlo  determined  values  of  M 


O 

n  »  3 

n  *  4 

n  -  5 

n  ■  10 

n.  20 

.01 

-.3138 

-.1948 

-.4976 

-.7315 

.1 

-.0500 

-.1171 

-.1732 

-.3324 

-.4630 

.3 

.0076 

.0143 

.0148 

.0363 

.0471 

.5 

.0708 

.1395 

.1924 

.3670 

.4830 

.7 

.1285 

.2516 

.3486 

.6478 

.85 69 

.9 

.1798 

.3487 

.4812 

.8821 

1.1807 

1.0 

.2030 

.3918 

.5392 

.9840 

1.3263 

1.5 

.2976 

•  5614 

.7630 

1.37  68 

1.9187 

2.0 

.3650 

.6760 

.9099 

1.6356 

2.3368 

5.0 

.5408 

•  9441 

1.2413 

2.2215 

3.3489 

10.0 

.6200 

1.0494 

1. 3684 

2.4426 

3.7  436 

- 

Table  12.  The  statistics  R,  R  and  R  as  functions  of  a  and  n 

- - - - - - * - - - - - — ■  o - — 


a 

n»5 

0 

H 

1 

« 

n  •  20 

0 

_  W 

H  R  R 

- 

ERR 

E  R 

R 

.0 

1.5123  1.4300  ^.3600 

1.3100  1.2780  1.2250 

1.1965  1.1800 

1.1630 

.1 

1.5428  1.4635  1.4070 

1.3379  1.3082  I.2550 

1.2200  1.2079 

1.1898 

.2 

1.5776  1.4995  1.4360 

1.3751  1.3454  1.2900 

1.2562  1.2461 

1.2252 

.3 

1.6166  1.5379  1.4670 

1.4215  1.3898  1.3300 

1.3051  1.2947 

1.2694 

.4 

1.6599  1.5787  1.500C 

1.4774  1.4414  1.3750 

1.3667  1.3537 

1.3224 

.5 

1.7075  1.6220  1.5350 

1.5426  1.5000  1.4250 

1.4411  1.4230 

1.3840 

.6 

1.7594  1.6677  1.5720 

1.6171  1.5658  1.4800 

1.5282  I.5027 

1.4644 

1.8156  1.7159  1.6110 

1.7009  1.6386  1.5400 

1.6280  1.5927 

1.5334 

I.876O  1.7665  1.6520 

1.7940  1.7186  1.6050 

1.7406  1.6931 

1.6212 

1.9408  I.8195  1.6950 

1.8965  1.8058  1.6750 

1.8658  1.8039 

1.7178 

2.0098  1.8750  1.7400 

2.0083  1.9000  1.7500 

2.0038  1.9250 

1.8230 

n  »  5 


R-  1.5123  +  0.2835  a  + 

0.2140  a2 

i 

a  «  -  0.6624+  2.162 

V  R  -  1.4184 

B -  1.4300  + 0.3230  0  + 

0.1220  a2 

1 

a  »  -  1.3238+  2.863 

V  B  -  1.2162 

R-  1.3800  +  0.2600  0  + 

0.1000  a2 

I 

a  -  -  1.3000+  3.162 

’/  R  -  1.2110 

n 

s 

10 

If- 1.3100  + 0.2321  a  + 

0.4662  a 

f 

a  -  -  0.2489  +  1.465 

f  R  -  1.2811 

R-  1.2780  +  0.2660  0  + 

0.3560  a 

> 

a  »  -  0.3736+  1.676 

V  H  -  1.2283 

B-  1.2250  +  0.27  50  a  + 

0.2500  a2 

1 

a  -  -  0.5500  +  2.000 

V  R  -  1.1494 

n 

s 

20^ 

B«  1.1965  +  0.1711  a+f 

0.6362  a 

; 

a  -  -  0.1345  +  1.254V  R  -  1.1850 

R  =  1.1800  +  0.2270  a  + 

0.5180  a2 

1 

a  -  -  0.2191  +  I.389V  B  -  1.1551 

R  *  1.1630  +  0.2240  a  + 

0.4360  a2 

i 

a  «  -  0.256" 

V  R  -  1.1342 
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Monte-Carlo  determined  valnes  of  n.VarfR^) _ 


0.3870 

0.4016 

0.4373 

0.4780 

0.5239 

0.5730 

0.5979 

0.7113 

0.7896 

0.7760 

0.5350 


n  ■  4 

n  »  5 

n  ■  10 

n- 20 

0.3881 

0.3555 

0.2112 

0.1510 

0.4197 

0.3926 

0.2430 

0.1748 

0.5061 

0.5012 

0.3456 

0.2512 

0.6209 

0.6650 

0.5470 

0.4196 

0.7716 

0.9107 

0.9741 

0.8742 

0.9580 

1.2538 

1.8275 

2.0646 

1.0617 

1.4601 

2.4847 

3.1782 

1.6115 

2.6895 

8.6676 

19.8124 

2.0745 

3.8639 

18.5044 

65.0248 

2.4740 

5.4802 

50.7443 

404.3828 

1.6978 

3.9759 

42.5263 

406.5855 

Table  1 


1.6350 

0.7443 

0.7424 

0.7376 

0.7280 

0.72J4 

0.6759 

0.6216 

0.3744 

0.2154 


1.3610 

1.3430 

1.3346 

1.3204 

1.2942 

1.2764 

1.1592 

1.0308 

0.5501 

0.2948 


5.8002 

1.8820 

1.8258 

1.8020 

1.7888 

1.7542 

1.7458 

1.6084 

1.4464 

0.6106 

0.4476 


3.9891 

3.5083 

2.9797 

3.0208 

3.3156 

3.6320 

3.7664 

4.1188 

4.0831 

2.6794 

1.5450 


6.6079 

4.9986 

3.5264 

4.0518 

5.4286 

7.1010 

7.9406 

11.3314 

13.0092 

10.6492 

6.3194 


Table  16 


+.  s 

.  The  variances  n.Var  a,  „  and  n.Var  #„,  as  functions  of  a  and  n 
. . .  ■ 1  1  ■  1  ■  id  ——————  ■■■  — . . . .  o  —  — 


n 

-5 

n 

-10 

n  - 

20 

n 

-  00 

*0 

n.Var  a+ 

n.Var  a* 

„  + 
n.Var# 

n.Var#* 

n . Var  #+ 

n.Var  a* 

n.Var  a+ 

n.Var#* 

.01 

4*2919 

5.3093 

3.6241 

0.9198 

4.4698 

0.8342 

_ 

0.8647 

.1 

3.6870 

1.8745 

2.2963 

0.8960 

1.9644 

0.6986 

- 

0.5507 

.3 

2.9539 

2.2243 

1.3194 

0.9748 

0.8220 

0.6308 

- 

0.4305 

.5 

2.6867 

2.7465 

1.1218 

1.3109 

0.6438 

0.9604 

0.1557 

0.8894 

mm 

2.6786 

3.5089 

1.2443 

1.9999 

0.7754 

1.7857 

0.3445 

2.0421 

EB 

2.8036 

4.6192 

1.5924 

3.2496 

1.1916 

3.4570 

0.7073 

5.1720 

\m 

2.8843 

5-3733 

1.8323 

4.2445 

1.5253 

4.8617 

1.0000 

8.0000 

Table  17.  The  efficiencies  Eff  a*„  and  Eff#f,  as  functions  of  a  and  n 
- . ""  . . . .  12  - - —  23 - — -  o - «— 


~  — 

“o 

— 

n  - 

5 

n  - 

10 

n  -  20 

n  - 

QD 

n.Var  & 

Eff  a+ 

Eff  a* 

Eff  a+ 

Eff  a* 

Eff  «+ 

Eff#* 

.0 

- 

— 

— 

- 

- 

- 

- 

53.30 

0.4765 

.1 

10.75 

21.14 

17.26 

44.23 

20.17 

56.72 

- 

72.79 

0.4008 

.2 

- 

- 

- 

- 

- 

- 

- 

88.34 

0.3304 

.3 

8.99 

11-94 

20.12 

27.24 

32.30 

42.09 

- 

61.63 

0.2653 

.4 

— 

- 

- 

— 

- 

- 

- 

34.63 

0.2058 

•  5 

5.66 

5.53 

13.55 

11.59 

23.60 

15.82 

97.59 

17.09 

0.1520 

.6 

- 

- 

- 

- 

- 

- 

93.19 

15.94 

0.2188 

.7 

11.12 

8.49 

23.94 

14.89 

38.41 

16.68 

86.46 

14.59 

0.2979 

.8 

- 

- 

- 

- 

- 

- 

78.32 

11.68 

0.3891 

.9 

17-56 

10.66 

30.92 

15.15 

41.32 

14.24 

69.62 

9.52 

0.4924 

1.0 

21.07 

11.31 

33.17 

14.32 

_ 

39.85 

12.50 

60.79 

7.60 

0.6079 
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